Understanding Pentaquark States in QCD 
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We estimate the mass of the pentaquark state with QCD sum rules and find that pentaquark states 
with isospin / = 0, 1, 2 lie close to each other around (1.55±0.15) GeV. The experimentally observed 
baryon resonance O"''(1540) with 5* = +1 can be consistently identified as a pentaquark state if its 
= i . Such a state is expected in QCD. If its parity is positive, this pentaquark state is really 
exotic. Now the outstanding issue is to determine its quantum numbers experimentally. 
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Quantum Chromodynamics (QCD) is believed to be 
the underlying theory of the strong interaction. In the 
high energy regime, QCD has been tested up to 1% level. 
In the low energy sector, QCD is highly nonperturbative 
due to the non-abelian SUc(3) color group structure. It 
is very difficult to calculate the whole hadron spectrum 
from first principles in QCD. With the rapid development 
of new ideas and computing power, lattice gauge theory 
may provide the final solution to the spectrum problem 
in the future. But now, people have just been able to 
understand the first orbital and radial excitations of the 
nucleon with lattice QCD in the baryon sector [1] . 

Under such a circumstance, various models which are 
QCD-based or incorporate some important properties of 
QCD were proposed to explain the hadron spectrum and 
other low-energy properties. Among them, it is fair to 
say that quark model has been the most successful one. 
It is widely used to classify hadrons and calculate their 
masses, static properties and low-energy reactions [2]. 
According to quark model, mesons are composed of a 
pair of quark and anti-quark while baryons are composed 
of three quarks. Both mesons and baryons are color sin- 
glets. Most of the experimentally observed hadrons can 
be easily accommodated in the quark model. Any state 
with the quark content other than qq or qqq is beyond 
quark model, which is termed as non-conventional or ex- 
otic. For example, it is hard for /o(980)/ao(980) to find 
a suitable position in quark model. Instead it could be a 
kaon molecule or four quark state [3]. 

However, besides conventional mesons and baryons, 
QCD itself does not exclude the existence of the non- 
conventional states such as glueballs {gg,ggg,- ■ ■), hy- 
brid mesons (qqg), and other multi-quark states {qqqq, 
QgggQj QQQQQQ^ iggggg-,' ' ')■ in fact, hybrid mesons can 
mix freely with conventional mesons in the large Nc limit 
[4] . In the early days of QCD, Jaffe proposed the H parti- 
cle [5] with MIT bag model, which was a six quark state. 
Unfortunately it was not found experimentally. 

In the past years there have accumulated some ex- 
perimental evidence of possible existence of glueballs 
and hybrid mesons with exotic quantum numbers like 
J-^*-^ = 1~+ [3]. Recently BES collaboration observed a 
possible signal of a proton anti-proton baryonium in the 
J/vt radiative decays [6]. But none of these states has 
been pinned down without controversy. 

Recently LEPS Collaboration at the SPring-8 facility 



in Japan observed a sharp resonance 8+ at 1.54 ± 0.01 
GeV with a width smaller than 25 MeV and a statisti- 
cal significance of 4.6cr in the reaction 7rt — > K'^K~n 
[7]. This resonance decays into K^n, hence carries 
strangeness S — 4-1. Later DIANA Collaboration at 
ITEP observed the same resonance at 1539 ± 2 MeV 
with a width less than 9 MeV in a different reaction 
K+Xe -> e+Xe' K°pXe' [8]. Now 6+ decays into 
K'^p. The convincing level is AAa. Very recently CLAS 
Collaboration in Hall B at JLAB also observed in the 
K^n invariant mass at 1542 ± 5 MeV in the exclusive 
measurement of the 'yd —>■ K^K~pn reaction [9]. The 
statistical significance is 5.3cr. The measured width is 
21 MeV, consistent with CLAS detector resolution. The 
isospin, parity and angular moment of the particle 
has not been determined rigorously yet. There is only 
preliminary hint that might be an iso-singlet from 
the featureless M{K'^p) spectrum in the CLAS measure- 
ment. 

Diakonov et al. proposed the possible existence of the 
S — 1 — resonance at 1530 MeV with a width 
less than 15 MeV using the chiral soliton model [10], 
which partly motivated the recent experimental search 
of this particle. The authors argued that 0+ is the light- 
est member of the anti-decuplet multiplet which is the 
third rotational state of the chiral soliton model. Assum- 
ing that the A^(1710) is a member of the anti-decuplet, 
0+ mass is fixed with the symmetry consideration of 
the model. Gao and Ma also discussed the experimental 
search of this pentaquark state [11]. Recently Polyakov 
and Rathke suggested possible photo-excitation of the 
pentaquark states [12]. 

Capstick, Page and Winston pointed out that identi- 
fying iV(1710) as a member of the anti-decuplet in the 
chiral soliton model is kind of arbitrary [13]. Instead, 
if the anti-decuplet Pu is A^(1440), the 9+ would be 
stable as the ground state octet with a very low mass 
while 0+ would be very broad with the anti-decuplet 
Pii being N{21QQ) [13]. Furthermore, if the decay width 
of the anti-decuplet A^(1710) was shifted upwards to be 
comparable with PDG values, the predicted width of 0"*" 
particle would have exceeded the present experimental 
upper bound [13]. The authors hypothesized that 6+ 
could be an isotensor resonance with negative parity [13]. 
The isospin violating decay width is naturally very small. 
Hence 9+ is a narrow state. Note similar isospin violat- 
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ing decay mechanism makes Dsj{2317) and £)sj(2460) 
states very narrow [14-16]. 

Stancu and Riska studied the stabiUty of the strange 
pentaquark state assuming a flavor-spin hyperfine inter- 
action between quarks in the constituent quark model 
[17]. They suggested that the lowest lying p-shcU pen- 
taquark state with positive parity could be stable against 
strong decays if the spin-spin interaction betwenn the 
strange antiquark and up/down quark was strong enough 
[17]. The group classification of strange pentaquarks is 
given in [18]. 

In Ref. [19] Hosaka emphasized the important role of 
the hedgehog pion in the strong interaction dynamics 
which drives the formation of the 0+ particle and the 
detailed energy level ordering. Hyodo, Hosaka and Oset 
suggested measuring the 6+ quantum numbers using the 
reaction if+p Tr+if+n [20]. 

In Ref. [21] Karlincr and Lipkin discussed the dynam- 
ics of the diquark-triquark pentaquark state with a rough 
estimate of the mass of 1592 MeV and I = 0, = 
in the constituent quark model. The same authors dis- 
cussed pentaquarks containing a heavy baryon [22]. 

Very recently Jaffe and Wilczek suggested that the ob- 
served 0+ state could be composed of an anti-strange 
quark and two highly correlated up and down quark pairs 
arising from strong color-spin correlation force [23]. The 
resulting of 0+ is | . 

Due to the low mass of 6+, we think its angular mo- 
mentum is likely to be one half. In this paper we shall 
employ QCD sum rules to estimate the mass of the pen- 
taquark states with J = ^ and 7 = 0, 1, 2. 

The method of QCD sum rules incorporates two ba- 
sic properties of QCD in the low energy domain: con- 
finement and approximate chiral symmetry and its spon- 
taneous breaking. One considers a correlation function 
of some specific interpolating currents with the proper 
quantum numbers and c;ak:ulates the correlator pertur- 
batively starting from high energy region. Then the res- 
onance region is approached where non-perturbative cor- 
rections in terms of various condensates gradually be- 
come important. Using the operator product expansion, 
the spectral density of the correlator at the quark gluon 
level can be obtained in QCD. On the other hand, the 
spectral density can be expressed in term of physical ob- 
servables like masses, decay constants, coupling constants 
etc at the hadron level. With the assumption of quark 
hadron duality these two spectral densities can be related 
to each other. In this way one can extract hadron masses 
etc. For the past decades QCD sum rules has proven 
to be a very powerful and successful non-perturbative 
method [24,25]. 

We use the following interpolating current for the 1 = 2 
pentaquark state where three quarks and the remaining 
qq pair are both in a color adjoint representation. We 
note in passing such a choice of color configuration is not 
unique [26]. 



r]2{x) = e°''"'[ul{x)C'y^Ub{x)]'y^'y5Ue{x)se{x)i'j5Uc{x) 

(1) 

for a charge Q = +3 state. Or 

mix) = -^{e'^'''[ul{x)Cj,Ukix)]ri5de{x) 

xseix)i-f5dcix) + {u^ d)} (2) 

for a charge Q = +1 state. The isospin of the current is 
shown in its lower index. a,b,c etc are the color indices. 
T denote transpose. C is the charge conjugation matrix. 
The property (C7^)^ = ensures the current is sym- 
metric under the exchange of the two quark fields inside 
the bracket. Both currents in Eq. (1) and (2) carry an 
isospin of two. 

We want to emphasize the isospin and color structure 
of these currents guarantee they will never couple to a 
K'^n molecule or any other K~^n intermediate states 
since the isospin of a K~^n system can only be zero or 
one. One may worry about KA intermediate states since 
their total isospin can be two. However the energy of KA 
intermediate states are greater than uik + tua = 1726, 
which is much larger than observed 0+ mass. In some 
cases, if the energy of the S-wave intermediate state is 
significantly lower than the resonance, its contribution 
could turn out to be quite important [27,16]. Angular 
momentum and parity conservation requires that only 
D-wavc KA intermediate states could contribute if 
of is ^ . When its is ^ , then the KA com- 
bination should be of p-wave. Either P-wave or D-wave 
KA continuum contribution is negligible compared with 
the lower 0+ narrow resonance contribution. So these 
currents couple to / = 2 strange pentaquark states. 

For 7 = pentaquark state the interpolating current 

is 

mix) = -^e''''-[ul{x)Cj5db{x)]{ue{x) 
v2 

xse{x)i'y5dc{x) - {uf-^ d)} (3) 

where ((775)"^ = — C'75 ensures the isospin of the up and 
down quark pair inside the first bracket to be zero. The 
anti-symmetrization in the second bracket ensures that 
the isospin of the other up and down quark pair is also 
zero. 

For 7=1 pentaquark state the interpolating currents 
are 

mix) = -^e''''[uJix)C-/^dbix)]{n5Ueix) 

XSe{x)ij5dc{x) - (m ^ d)} (4) 

or 

mix) = -^e«^^[<(a;)C754W]{«eW 

xseix)i'y5dcix) + iu f-^ d)} (5) 
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In the following we use currents (2), (3) and (5) to per- 
form the calculation. 

We introduce the overlapping amplitude fj of the in- 
terpolating currents with the corresponding pentaquark 
states. 



|ry,(0)b,/ = z)=/,ii(p) 



(6) 



where u{p) is the Dirac spinor of pentaquark field with 
I = J- 

We consider the correlator 



\T{7j,{x),fj,io))\o) =n(p)p + n'(p) (7) 



where fj = 77^70 and p = -7^. Throughout this note we 
focus on the chirality even structure n(p). At the hadron 
level it can be written as 



n(p) 



■'3 



Mf 



higher states 



(8) 



where Mj is the pentaquark mass. 

On the other hand, it will be calculated in terms of 
quarks and gluons. For example, for the 1 — 2 case, the 
correlator reads 



xTr [175 • iSl>^{-x) ■ ij5 ■ iScc'{x)] j'^js ■ iSee'{x) ■ 7'' 75 



-2e 



abc b' c 



Tr [iSaa'ix) ■ 7i.C • iS'[f,,{x) ■ C"f^ 



X7^75 • iSec'ix) ■ «75 • iS^,^{-x) ■ 175 • iSce'{x) ■ Yl5 



+4e-b-e-'f'c'^M^5 . iS,a'{x) ■ iM^Sbc'ix) ■ ^75 x 
iSl,^{-x) ■ ij5 ■ iScb'ix)]'^Cjf^ ■ iSae'ix) ■ J^js} 

where iS^i^{—x) is the strange quark propagator in the 
coordinate space. 

After making Fourier transformation to the above 
equation and invoking Borel transformation to Eq. (7) 
we get 



e T- 



'Pjis) 



(9) 



where is the strange quark mass, Pj{s) is the spectral 
density and sq is the threshold parameter used to sub- 
tract the higher state contribution with the help of quark- 
hadron duality assumption. Roughly speaking is 
around the first radial excitation mass. The spectral den- 
sity reads 



P2 



46717^8 2i36!7r' 



15367r4^3 



Pi 



5s^ 



,13 



2i77!7r8 15367r4 48 



1 



(qq) +o{qq){ss)] 



3s5 



Po 



4»71tt» 15367r4 42 



{qq)' 



11 

24 



{qq){ss)] 



U32 



{qqy- 



ss 



864 



{qq)'m 



where we have used the factorization approximation for 
the multi-quark condensates. 

In order to extract Mj we first take derivative of Eq. 
(9) with respect to l/T^. Then we divide it by Eq. (9) 
to get 



Mf 



(10) 



with p'{s) — sp{s) except that p'{s) does not contain the 
last term in p{s). 

In the numerical analysis of the sum rules the values 
of various QCD condensates are (ss) = —(0.8 ± 0.1) * 
(0.24 GeV)3, (g^GG) = 0.48 GeV. We use m,(lGeV) = 
0.15 GeV for the strange quark mass in the MS scheme. 

Numerically we have M2 = (1.53 ± 0.15)GeV, Mi = 
(1.59 ± 0.15)GeV, Mq = (1.56 ± 0.15)GeV, where the 
central value corresponds to T = 2 GeV and sq = 4 
GeV^. The variation of M2,Mq with both T and sq is 
shown in Figures (l)-(2), which contributes to the er- 
rors of the extracted value, together with the truncation 
of the operator product expansion, the uncertainty of 
vacuum condensate values and the factorization approx- 
imation of the multiquark condensates. In the working 
interval of the Borel parameter Mj is reasonably stable 
with T. However the level ordering of the qentaquarks 
with isospin should not be taken too seriously due to the 
large uncertainty. 
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FIG. 1. The variation of the / = 2 pentaquark state mass 
M2 with the Borel parameter T (in unit of GeV) and the 
continuum threshold so- From bottom to top the curves cor- 
respond to So — 3.61,4.0,4.41 GeV^ respectively. 

Our results suggest that strange pentaquark states 
with / = 0, 1, 2 lie close to each other, which is very dif- 
ferent from the prediction of chiral soliton model where 
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only the lightest 7 = member is around 1540 MeV and 
all other states are several hundred MeV higher [10]. As 
we mentioned before, our interpolating currents can cou- 
ple to pentaquark states with either positive or negative 
parity. QCD sum rules approach only picks out the low- 
est mass pentaquark state in the specific channel. It is 
impossible to determine its parity in our formalism. 
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FIG. 2. The variation of the / = pentaquark state mass 
Mo with the Borel parameter T (in unit of GeV) and the 
continuum threshold so- From bottom to top the curves cor- 
respond to So — 3.61,4.0,4.41 GeV^ respectively. 

One may tend to think that O is a = i strange 
pentaquark state from previous quark model experience. 
One radial or orbital excitation of constituent quarks typ- 
ically carries energy around 400 — 500 MeV. Naively one 
would expect the strange pentaquark mass with positive 
parity to be around 4Af„ -I- + 400 2050 MeV, where 
Mu = 300 MeV and Ms = 450 MeV are roughly the 
up and strange quark constituent mass. The additional 
~ 400 MeV comes from one orbital excitation for the 
positive-parity If future experiments further estab- 
lish the O"*" parity to be positive, then this particle is 
really very exotic. Otherwise, a = ^ O"*" particle 
is an expected strange pentaquark state in QCD. The 
outstanding issue is to measure the quantum numbers of 
the O^ particle, especially the parity, angular momen- 
tum and isospin in the future experiments. 

As a byproduct, we want to mention that the radial 
excitation 8' of the strange pentaquark may lie around 
2.0 GeV. If 8 is an iso-tensor state as suggested in 
[13], 8' should mainly decay into final states contain- 
ing one or two pions like 8' — > K^A n'^K^n, 8' — > 
87r7r — > K^mrn, which conserve isospin symmetry. The 
isospin conserving decay mode 8 — > n'^K^n is kine- 
matically forbidden for an isotensor pentaquark. If its 
isospin is or 1, then the dominant mode should be 
8' K^n, 8' Qtttt — > K^rnnr. In other words, 
the decay mode 8' n'^K^n with a single pion in the 
final state can be used to distinguish the isospin of the 
strange pentaquark if enough events are collected in the 
future experiments. 
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